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I. INTRODUCTION 

The random-field Ising Magnet (RFIM) is a prototyp- 
ical model for magnetic systems with quenched disorder. 
For d = 3 and higher dimensions,!^ it is known to undergo 
a second-order phase transition^ini at a critical tempera- 
ture Tc or disorder strength he {Tc{h) 4^ hc{T)): For low 
temperatures and weak disorder the ferromagnetic inter- 
actions dominate and the system is long-range ordered. 
For large temperature or strong disorder, the RFIM ex- 
hibits no long-range order and behaves like a paramagnet 
in a field. 

The quenched disorder used in earlier studies of the 
RFIM was mostly uncorrelated (^-correlated) .'^^ This 
is quite common in the literature when studying dis- 
ordered systems like percolation, random ferromagnets, 
spin glasses or polymers in random media. Nevertheless, 
real systems are always emerging from physical processes, 
hence correlations are present, which could play an im- 
portant role for its behavior. Here, we consider a tun- 
able, scale- free (power-law), i.e. long range, correlation 
to the random field to explore its influence on the crit- 
ical behavior. Please note that for an exponentially de- 
creasing correlation strength with a typical length scale 
S, via renormalizing the system beyond S, the behav- 
ior of the uncorrelated system should be recovered. The 
0{n) random- field model with long-range correlated dis- 
order was studied recently^* via functional renormaliza- 
tion group methods around d — A and for values n > 3, 
i.e., without including the Ising case n = 1. For other 
types of random systems, there exist already some stud- 
ies for the case of long-range correlated disorder^ e.g., for 
percolation,'^^ the diluted Ising ferromagnet,^^ random 
walkSj'^or elastic systems.'^ 

Now, we state our model in detail. The RFIM consists 
of Ising N = spins 5*^ = ±1 located on the sites of 



a cubic lattice with periodic boundary conditions in all 
directions. The spins couple to each other and to local 
net fields. Its Hamiltonian reads 

H^-JY, S^SJ - (% + H)S,. (1) 

It has two contributions. The first covers the spin-spin 
interaction, where J is the ferromagnetic coupling con- 
stant between two adjacent spins and («, j) denotes pairs 
of next-neighboured spins. The second part of the Hamil- 
tonian describes the coupling to local, and global fields 
h-qi and H, respectively. The factor h is the disorder 
strength used to trigger the phase transition. The global 
field is included only for technical reasons to calculate 
the susceptibility in the limit _ff — 0. The quenched lo- 
cal fields "qi are Gaussian distributed with zero mean and 
unity width. The important property of these fields is 
their spatial long-range correlation. It decays as a power 
law 

^(^^ ^ (^E^(^)^(^ + ^1) - (2) 

X 

with a tunable, well defined decay exponent a. The sym- 
bol (...) denotes the average over the quenched disorder 
and X are the positions of a lattice sites i. 

We will study in particular the values a — 
{— 1, — 2, — 3, — 7}. First of all, it is interesting to know 
whether this type of disorder is relevant with respect to 
the ordered case. A hint to the answer of this ques- 
tion comes from the case of systems with "random- 
temperature disorder", like the diluted ferromagnet: For 
a d-dimensional system, if \a\ > d, i.e. when the disor- 
der correlation vanishes rather quickly, the usually Harris 
criterion applies.'^ The Harris criteriorPSl states that the 
disorder is relevant \i dv — 2 < Q, v being the critical 
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exponent of the ordered system. For the d = 3 ferro- 
magnet, we have v = 0.6294(5) from Ref. ETJ hence the 
disorder is relevant, as known from the case of uncorre- 
lated disorder. For |a| < d, in particular for a — —\ and 
a = —2 as studied here, the disorder is relevant accord- 
ing to Ref. Hi for 2/|a| > v. Since 2/1 = 2 > 0.6294 
and 2/2 = 1 > 0.694, the disorder will be relevant also 
for these values of a. The results we present below in 
this work, although for a different type of disorder, are 
compatible with these predictions. Furthermore, we will 
consider the question whether the correlated disorder is 
different from the behavior of the uncorrelated disorder 
case. Our results show that the most exponents are com- 
patible within error bars with the values of the standard 
RFIM, but the combination ^jv shows a clear signature 
of non-universality. This is similar to the diluted Ising 
model, where the long-range correlated dilution clearly 
changes some but not all critical exponent^^ with re- 
spect to the uncorrelated cas,^^ 

The paper is organized as follows: In section |ll] we 
sketch the idea how to calculate Gaussian distributed cor- 
related random numbers. After that a brief description of 
the numerical ground-state approach is given. The mea- 
sured quantities and the methods to analyze the data are 



displayed in section III Our numerical results are pre- 
sented in section HVl Based on the results we discuss the 
extremes of correlated disorder. The last section contains 
the discussion and conclusion. An extensive summary of 
this work (1/10 of the length) can be downloaded from 
the Papercore database.'^^. 



II. NUMERICAL METHODS 

In this section, we first explain how we generated the 
samples of the correlated disorder. Second, we briefly 
outline the numerical approach to calculate the exact 
ground states of these samples. 

To obtain a realization of correlated random fields, we 
basically apply the ideas of Refs. I15|16|19l The recipe is, 
to demand for a convolution kernel <I>(r) which convolves 
iid random numbers u(r), such that 7j{r) = $(r) * u(r^ — 
^ix)u{f— x) (using periodic boundary conditions for 
u{.)) show a desired two point correlation. Power law 
correlations are created, using 



C(r)= (l+H2)a/2 



(3) 



The long range behaviour is the same as of a pure power 
law without a si ngular ity at the origin. This avoids zero- 
mode divergence .122251 

In Fourier space, the transformation J- given through 

fj{k) = '^ge^'^'^v{^)i the correlation function is equiva- 
lent to the spectral density. Applying the definition of 
C(r) from Eq. ^ results to 
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A convolution in real space turns to a multiplication in 
Fourier space: 



fj{k) = ^{k)u{k). 



(5) 



Now insert Eq.([5| into Eq.Q to determine the convolu- 
tion kernel. 

NC{k) = (^{k)u{k)^*{k)u*{k)) = |l>(fc)|2(|u(fc)|2). (6) 

We choose the real space random numbers u{x) as being 
distributed iid according a Gaussian with zero mean and 
variance one, 



{u{x)) = 
{u{x)u{y)) = 53,y, 



(7) 
(8) 



such that the variance (|{i(fc)p) N. This results in 
C{k) — |$(fc)p, such that we can choose the correlated 
random numbers (in the Fourier space). 



(9) 



The back transformed correlated random numbers ?7(r) 
are real numbers. Since u{r) G M, u{—k) = u*{k). From 
Eq. Q we infer C{k) = C{\k\) e M+. So the back trans- 
formation Tjir) = {fj{k)){r) E M. 

For maximum flexibility, to test different types of 
correlation functions, we have implemented the Fourier 
transformation numerically, using the Fastest Fourier 
Transform in the West (FFTW) library, version 
An example of a realisation of the correlated disorder is 
shown in Fig. [ijfor different values of a. 



a=-7.0 




--NCik) . 



(4) 



FIG. 1: (color online) Slices of the same correlated disorder of 
a 97'^ lattice for different correlation strengths. The random 
fields are heat coded. Bright means high positive field and 
dark means high negative field. 



We tested our procedure by generating 20 realizations 
for a = —1 of the correlated disorder, calculating the 
two-point correlation Eq. ^ directly, and fitting Eq. ^ 
with variable exponent a. The correlation, shown in Fig. 
[2] and the resulting value a — 1.003(2) show that except 
for very small correlation at large distance the procedure 
works very well. 

Next, we mention shortly how the exact ground states 
are calculated. The phase space of the uncorrelated 
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FIG. 2: (color online) Mean of the two-point correlations for 
20 arbitrary samples of long-range correlated random fields 
for L = 49, a = — 1 (black line) and error (gray background). 
The dashed (magenta) line is a fit according to Eq. (|3| with 
a = 1.003(2) . The inset shows the same data as log-log-plot. 



since it is possible to calculate exact ground states at 
T = in a very efficient way for system sizes as large 
as = 141^ spins. Within this approactPsEH^ each re- 
alization of the correlated net fields {{rji],H) has to be 
mapped to a graph with -f 2 nodes and 2N + 1 edges 
with suitable edge capacities. On this graph a sophis- 
ticated maxi mum flow/minimum cut algorithm can be 
applied. 1^21311 The resulting minimum cut directly corre- 
spond to the GS spin configuration {Si} of that specific 
realization of the net disorder. We used the efficient 
maximum-flow subroutines implemented in the LEDA 
library.l^Sl 



III. QUANTITIES OF INTEREST 

From a GS spin configuration, some quantities of in- 
terest can be obtained directly, as the magnetization per 
spin 



1 



TV 



(10) 



RFIM consists of a ferromagnetic and a paramagnetic 
phase, see Fig. |3] The transition from one phase to the 
other can be triggered by varying the disorder strength 
h or the temperature T. Changing both along a path 
f{h,T) in the phase space leads to a critical point Pc — 
{hc^Tc) f{^h.T)- From renormalization group calculations, 
the RFIM is known^^ to exhibit the same critical be- 
havior at any Pc, except the temperature-driven phase 
transition point of the standard non-random Ising model. 
Hence, it is possible to focus on T = = const, and vary 
just h, to study the critical behavior along the full tran- 
sition line. We do not know a-priori, whether the phase 
diagram for the correlated-disorder case has the same 
property, nevertheless, it makes sense to concentrate, at 
least for our study presented here, also on T = 0. From 




FIG. 3; (color online) Schematic diagram of the phase space 
of the RFIM. The path f(h, T) (red line) shows an arbitrary 
path to cross the phase boundary. The small arrows denote 
the renormalization group flow. 

the computational point of view this is very favorable. 



and the bond energy per spin 



(11) 



Using these individual values, we calculate averaged 
quantities like the average magnetization m = (M). This 
disorder average (...) is performed always for a fixed 
value of h. We also consider the Binder cumulantP^ 



g{h,L) 




(12) 



A specific-heat-like quantity C{h) can be calculated as 
the numerical derivative of Ej with respect to h (see 
Ref. [5] for details). From here on we will refer to it as 
specific heat. 



C{h) 



{dEjjh)) 
dh 



(13) 



We also calculated the zero-temperature susceptibility 
dmih,H) 



dH 



(14) 



H=0 



as linear response of the magnetization to small homo- 
geneous magnetic fields H. Therefore, we apply small 
homogeneous fields at equidistant values Hi , 2_ffi , iHi 
and fit parabolas as function of H to the magnetizations 
m{h, H = 0), m{h, Hi), m{h, 2Hi) and m{h, 3Hi). For a 
fixed value of the disorder strength h the linear coefficient 
corresponds to the susceptibility x(^)- 

We will see that the results are compatible with second 
order phase transitions, such that the measured quanti- 
ties show power-law behavior close to the phase transition 
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point. To determine the critical exponents, we use the 
standard scaling forms, i.e., 



g{h,L) 




9 






(15) 


m(h, L) 




m 


{{h- 




(16) 


x{h,L) 




X 






(17) 


C{h,L) 




c 


((/^- 




(18) 



and apply a finite-size scaling analysis. For the Binder 
cumulant and the magnetization we use a nice tool which 
performs data collapses automatically;^^ It is based on a 
simplex algorithm and is written in python. 

The specific heat and the susceptibility show a maxi- 
mum close to the critical point, at some argument / of 
the universal functions x(0 and C{-)- Note that the peak 
positions for specific heat and susceptibility of the same 
system sizes L usually differ. Thus, also the value of / 
(and even the sign) may differ. From Eqs. (17) and (18 1 



it follows that the finite-size dependence of the positions 
of the maxima, respectively, scale as 



h*{L) = h,, + f -L-^'". 



(19) 



Furthermore, right at h*{L), the height of the maxima 
should scale as L^'" and 1°''" . In the case of a = other 
forms like a logarithmic divergence or a convergence to a 
constant ("cusp ") ha ve been observed for other systems 
in the literature.'^'^ Below we present the results we ob- 
tained for the position and the height of the peaks and 
test their scaling behaviour according to these scaling as- 
sumptions. 

As mentioned above, these quantities are average val- 
ues. They are strongly dependent on the set of disorder 
realizations taken into account. Hence, we perform an 
average usually over many thousands of realizations. We 
estimate the variability of these average values from 200 
bootstrap samples"^^ and quote this as error. 



IV. RESULTS 

We performed exact ground-state calculations for 
three-dimensional RFIMs for correlation strengths a — 
{— 1, — 2, — 3, — 7}. We considered system sizes ranging 
from L = 7 to L — 141. The number of disorder real- 
izations per system size and correlation strength can be 
found in Tab. [ij The actual number of calculated ground 
states is four times larger, since 4 different external fields 
are needed to obtain a susceptibility. The values of Hi 
are stated in the very right column of Tab. |lj 

Since the Binder cumulant exhibits no clear crossing 
(see Fig. |4| one could suspect that no phase transition 
is present. This is not the case as well will see in the 
following. To determine the phase transition points, we 
start by considering the average specific heat. In Fig. [5] 
the results can be seen for a = —7 to a = —1. As for 
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TABLE I: Number of disorder realizations per system size and 
correlation strength in thousands (10^) and the external field 
Hi, used to calculate the susceptibility. 
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FIG. 4: (color online) The Binder cumulant for a 
different system sizes. 



-3 and 



the uncorrelated RFIM, peaks can be observed clearly, 
which give evidence for the existence of a phase transi- 
tion also for the correlated case. We estimate the peaks 
by fitting parabolas over different intervals close to the 
maximum (for every bootstrap sample). The positions of 
the peaks in Fig. [5] move from right to left for increasing 
system sizes. To obtain the infinite-size limiting value he 
and an estimate for the critical exponent v of the cor- 
relation length, we fit the positions of the peaks to Eq. 
( 19 1, resulting in fit values as shown in the upper part of 
Tab. [n] Note that when determining the error bars from 
model fitting, we usually have not only taken the statis- 
tical error obtained from the fit routine (of the gnuplot 
program) but we have always also varied the range of 
sizes, to get an impression of possible systematic errors. 

In Fig. [5] for the case a — ~1, the peaks move very 
close to h = and the result from the fit for he is also 
close to zero. Therefore, another sensible ansatz is to set 
he — 0. Fit parameters for this ansatz are also shown 
in Tab. In] in the lower part. Both models are plotted 
in Fig. |6] as solid or broken lines, respectively. For a < 



5 




FIG. 5: (color online) Specific-heat-like quantity for different 
correlation strengths a = —7... a = —1 and system sizes 
L = 7, 11, 15, 21, 25, 35, 49, 69, 97, 117 and 141. For all values 
of a the curves appear in monotonic ordering, meaning the 
curve of L = 7 is on the very right and L = 141 is on the very 
left, as labeled on the upper left subplot. The lines are guides 
to the eyes only. 



— 1, clearly /ic > is better compatible with the data 
points, while for a = — 1 no real decision can be made 
here between /i^ > and = 0. Nevertheless, hc{L) 
also looks slightly curved in the double-logarithmic plot, 
hence he > appears more likely here as well. 





a = -7 


a = —3 


a = -2 


a = ~l 


he 


1.962(4) 


0.992(8) 


0.537(8) 


0.183(9) 


f 


3.2(1) 


6.0(3) 


7.5(2) 


9.8(3) 




0.79(2) 


0.84(2) 


0.81(2) 


0.84(2) 


assuming no critical point <4> /ic = 


h 


2.62(5) 


2.61(8) 


3.1(1) 


5.5(3) 




0.055(5) 


0.187(8) 


0.32(1) 


0.58(1) 



TABLE II: Fit parameters of Eq. [19] for the peak position of 
the specific-heat-like quantity. The upper part contains the 
parameters for finite h^. For the lower part /ic = was fixed. 



To determine the critical exponent a according to Eq. 



( 18 1, we analyzed the peak heights of the specific heat as 
shown in Fig. [7j They increase up to L « 50 for all cor- 
relation strengths a and decrease for larger L. Thus, no 
clear scaling is visible. This could be due to very strong 
finite-size corrections. Therefore, under the assumption 
that the specific heat decreases in a power-law fashion, 
we fitted the data points for very large system sizes a 
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FIG. 6: (color online) Peak positions of the specific-heat-like 
quantity as function of the system sizes for a — ~7 . . .a = 
— 1. The solid lines are fits assuming a finite critical value 
according to Eq. |19[ the broken ones are fits imply he — 0. 



power law of the form 

C{h,L) ^ kL"/" . 



(20) 



The achieved exponents are small and negative can be 
found in Tab. Ill On the other hand it may be that 
the specific heat levels off for even larger system sizes, 
which would give the leading behavior a = 0. In Sec. |VI[ 
we will discuss these two options in connection with the 
Rushbrooke inequalitjBSl and see that a = appears to 
be more likely. 





a = ~7 


a = —3 


a = -2 


a = -1 


k 


2.3(2) 


1.7(1) 


1.9(1) 


1.0(1) 


ajv 


-0.04(2) 


-0.11(2) 


-0.26(1) 


-0.24(3) 



TABLE III: Fit parameters of Eq. [20] of the height of the 
maxima of the specific-heat as shown in Fig. [7] 

We now turn to the susceptibility. The phase tran- 
sition is signaled by a divergence of the susceptibility. 
An increasing peak can be seen for a = — 2 in Fig. [8] as 
example. The peaks are estimated in the same way as 
for the specific heat. The resulting maxima are tuples 
(h* (L) , XmaxiL)) of position and height. 

For the peak position of the susceptibility we assumed 
the same model as we did for the specific heat. The mod- 
els and data points can be found in Fig. [9]ln particular 
for a — —1, the error bars are quite large, despite the 
large number of samples, which is for the largest system 
sizes considerably higher compared to the cases a < —1. 
To understand this behavior we studied the degree of 
non-self averagin^^II and we calculated 



R^iL) = var(x(L))/(x(i))' 



(21) 



We found to stay approximately constant for increas- 



ing L, as shown in Fig. 10 for susceptibility measured at 
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FIG. 7: (color online) Peak heights of the specific heat as 
function of the system sizes for a = —7 ... a — —1. The lines 
are fits assuming a power law decay for large system sizes. 



FIG. 9: (color online) Peak positions of the susceptibility as 
function of the system sizes for a = —7 . . . a = —1. The solid 
lines are fits assuming a finite critical value according to Eq. 
[191 the broken ones are fits imply he = 0. . 
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FIG. 8: (color online) Susceptibility for a — —2. 



rule out a he = for a = — 1, 
clearly h^, > holds. 



while for a < — 1 again 



the peak positions. We found the same behavior qualita- 
tively for different fixed values of h, which shows that the 
correlated RFIM is non-self averaging for a large range of 
the disorder parameter, as many other systems exhibiting 
quenched disorder. In particular, the results in Fig. |10| 
show that the degree of non-self averaging is strongest for 
a = — 1, which explains the large error bars. To achieve 
much smaller error bars for the susceptibility, a much 
larger number of samples would be necessary, which is 
beyond the capacity of our numerical resources. 

For the finite-size scaling of the peak-positions, we 
tested Eq. ( 19 1 using the saturating ansatz {h^ included 
in the fit) as well as a pure power law decay (via hc = 0). 
The fit parameters for both models can be found in Tab. 
|IV[ Again the saturating model brings up, within the 
present accuracy, the same infinite-size critical point he 
as we found before. Due to the error bars, we can not 





a = -7 


a = —3 


a = -2 


a = -1 


he 


1.97(1) 


0.98(4) 


0.56(4) 


0.20(5) 
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4.8(5) 


7(1) 


9(1) 


11(2) 


l/v 


0.83(5) 


0.84(8) 


0.85(7) 


0.88(8) 


assuming no critical point <4> /ic = 


/o 


3.06(8) 


3.4(3) 


4.5(4) 


7.7(7) 




0.088(7) 


0.24(2) 


0.41(3) 


0.66(3) 



TABLE IV: Fit parameters of Eq. [T9]of the peak position of 
the susceptibility. The upper part contains the parameters 
for finite he. For the lower part he = Q was fixed. 

In contrast to the specific heat, the peak heights of the 
susceptibility shows a clear power law behaviour for all 
studied correlation strengths, see Fig. 
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Thus, in the 

thermodynamic limit the susceptibility diverges. Com- 
pared to the peak positions displayed in Fig. [9j the fluc- 
tuations for the peak height are much smaller, thus a 
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a = -7 


a = —3 


a = -2 


a = -1 


bo 


0.064(3) 


0.049(2) 


0.049(2) 


0.041(3) 




1.56(1) 


1.45(1) 


1.34(2) 


1.20(2) 



TABLE V: Fit parameters of the peak heights of the suscep- 
tibiUty when fitting according to boL^'''^ . 
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FIG. 11: (color online) Peak heights of the susceptibility as 
function of the system sizes for o = — 3...a = — 1. 



clear power-law behavior is visible. The critical expo- 
nents 7/1^, as obtained from a power- law fit, are displayed 
in Tab. [V] The values decreases with increasing a. 

For a finite-size analysis of the Binder cumulant and of 
the magnetization we performed data collapses according 
to Eqs. (15) and (16). Example data collapses for the 



magnetization and as inset for the Binder cumulant for 
a 



-3 are shown in Fig. 12 One sees that the quality of 
the collapses is very good. These lead to sets of critical 
values and exponents as shown in Tab. IVl] 
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FIG. 12: (color online) Data collapse of the magnetization 
and of the Binder cumulant (inset) for a = — 3. 





a = -7 


a — —3 


a = -2 


a = -1 


he 


1.94(1) 


0.95(2) 


0.47(4) 


0.12(3) 


l/u 


0.78(2) 


0.77(4) 


0.75(3) 


0.78(4) 


fijv 


0.005(5) 


0.03(2) 


0.01(1) 


0.01(8) 



TABLE VI: Critical value he and correlations length exponent 
v derived from the finite size scaling analysis of the Binder 
cumulant and the magnetization. These values are obtained 
via data collapses. 



V. MINIMUM AND MAXIMUM 
CORRELATED DISORDER 

The minimum correlation range of the disorder is a — >■ 
—00 = (5-correlated disorder, i.e. the normal RFIM. The 
other extreme is a — >■ 0~ . As it is illustrated in Fig. [l] 
by increasing the correlation strength a, the regions of 
sites with almost the same sign of the field get larger and 
larger, while keeping the same average close to 0. Hence, 
for a — >■ 0~ one can imagine each realization of the dis- 
order being bi-parted. Bi-parted means, to find two dis- 
tinct clusters A = {r]i > 0} and B — {r]i < 0}, see Fig. 
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An Imry-Ma type of argument would read as follows: 
In such a disorder realization, for a state where all spins 
are aligned with its local field, the resulting interface en- 
ergy between the clusters A, B would be Ei ^ L'^f with 
fractal exponent d — 1 < df < d. This competes against 
the field energy Eh ~ —hL'^: If Ej < \Eh\, the ground- 
state will be an ordered phase, otherwise both clusters are 
locally aligned. A T = phase transition occurs when 
\Efi\ = Ej. From this we see that the finite-size critical 
point scales in the limit a ^ 0~ a.s he ^ L'^f-'^ — > 
for L ^ 00. This means for highly correlated but ar- 
bitrarily small disorder the RFIM will behave as a su- 
per paramagnet in a field in the thermodynamic limit. 
Though, there would be no disorder-driven phase tran- 
sition anymore. Furthermore, since the number of spins 
at distance r scales as r'^~^ but the disorder correlation 
decreases only as r", its appears plausible that even for 
a finite range of a < values the cumulative effect of the 
correlation might dominate and indeed he = already 
for these values of a. This explains why the result for 
a = — 1 is ambiguous. 



VI. CONCLUSION AND DISCUSSION 

We have presented the results of exact ground-state 
calculations of the RFIM with correlated disorder for 
different correlations strengths. To numerically calcu- 
late the ground states, we have applied a mapping to 
the maximum-flow problem. Using efhcient polynomial- 
time-running maximum- flow/minimum-cut algorithms, 
we were able to study large systems sizes up to = 141'^. 

We studied different quantities like magnetization. 
Binder cumulant, susceptibility and a specific heat-like 
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quantity and applied finite-size scaling techniques to ob- 
tain the critical exponents. The combined results for 
the critical exponents are shown in Tab. |VII| We tested 
the two possibilities for the values of a, by applying the 
Rushbrooke inequality a + 2/3 + j > 2 which holds usu- 
ally as equality.!^ When choosing a = 0, the Rushbrooke 
equation is fulfilled in all cases within error bars. For 
the values of a quoted in Tab. |III[ obtained via fitting 
the data for just the few largest systems sizes, the Rush- 
brooke sum is (assuming he > 0) considerably smaller 
than 2 for a — —2,-1. Hence, the value a — appears 
to be more likely. Note in all cases, the values quoted 
in the table are compatible within error bars with the 
results for the uncorrelated case, in particular due to the 
relative large error bar for the critical exponent 7. Never- 
theless, the data for the peak heights of the susceptibility 
(Fig. 11 ) show a trend towards a smaller slope when in- 
creasing a from —7 to —1: The results for ^/v, which 
are very precise (see Tab. [v| are clearly different within 
error bars. In this case, to still fulfill the Rushbrooke 
inequality, the true value for in particular for values 
a — —2 and a = — 1, should be larger, at or somehow 
above the upper bounds given the standard error bars. 
Hence, it is quite likely that the correlation of the dis- 
order creates non-universality for the RFIM, as in the 
case of the diluted ferromagnet.'^^. Also, among our re- 
sults, the case a = — 1 is special. As discussed above, it 
appears plausible for a = — 1 that the cumulative effect 
of the correlation might dominate. This would lead to 
he = already for a = —1. The result /i^ = is in- 
deed possible, as shown in last row of Tab. |Vn| Also on 
the base of the Rushbrooke sum, we cannot take decision 
on this issue. Nevertheless, numerically the inequality is 
much better fulfilled when assuming h,, ~ 0. 



(Grofirechner OLdenburg fiir Explizit Multidisziplinare 
Forschung) and the HERO (High-End Computing 
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FIG. 13: (color online) System with bi-parted random field 
(with periodic boundaries). 
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V 


RS sum 


—00 


2.27(1) 


~ 


0.016(7) 


2.1(1) 


1.37(9) 


2.1(2) 


-7 


1.96(2) 





0.01(1) 


2.0(3) 


1.26(8) 


2.0(3) 


-3 


0.97(2) 





0.04(6) 


1.8(4) 


1.2(2) 


1.8(5) 


-2 


0.52(5) 





0.01(3) 


1.7(3) 


1.2(2) 


1.7(4) 


-1 


0.17(5) 





0.05(10) 


1.5(3) 


1.2(2) 


1.5(5) 


-1* 








0.1(1) 


1.9(1) 


1.61(9) 


2.0(4) 



TABLE VII: Final results for the value of ftc, the critical 
exponents a, /3, 7 and v and the Rushbrooke (RS) sum a -I- 
2/3 -I- 7. The first fine shows the result for the uncorrelated 
case taken from Refs. i6i7i The last line is for the assumption 
/ic(a = -1) = 0. 
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